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Abstract 


A numerical rate equation model for a continuous wave iodine laser 
with longitudinally flowing gaseous lasant is validated by approximating 
two experiments that compare the perfluoroalkyl iodine lasants n-CjFyl 
and t-CjFpI. The salient feature of the simulations is that the production 
rate of the dimer ( C 4 F 9 ) 2 is reduced by one order of magnitude relative to 
the dimer (CjFy) 2 . The model is then used to investigate the kinetic effects 
of this reduced dimer production — especially how it improves output 
power. Related parametric and scaling studies are also presented . When 
dimer production is reduced, more monomer radicals (t-C 4 F 9 ) are avail- 
able to combine with iodine ions , thus enhancing depletion of the laser 
lower level and reducing buildup of the principal quencher, molecular 
iodine . Fewer iodine molecules result in fewer downward transitions from 
quenching and more transitions from stimulated emission of lasing pho- 
tons. Enhanced depletion of the lower level reduces the absorption of las- 
ing photons. The combined result is more lasing photons and 
proportionally increased output power . 
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lasant heating from pump, W * m 
quenching rate coefficient, cm 3 • sec 
gas constant, 8.314 J ■ mol -1 • K“ l 


R 


Mx 10“ 


J • kg - ' ■ K 1 


radical n-C 3 F 7 , t-C 4 F 9 , or f-C 3 F 7 
reflectivity of highly reflecting mirror, 3 .0 
reflectivity of output mirror, 0.75 
radial coordinate in cylindrical frame, m 

f 

normalized radial coordinate, — 

inner radius of laser tube, m 

distance along inner perimeter of laser tube, m 

arc length, m 

_2 

solar constant, 1 .35 kW • m 
lasant temperature, K 
lasant flow speed 

axial distance downstream of pump entrance 


v p pump frequency, Hz 

^ photodissociation rate of [RI], sec -1 

photodissociation rate of [I 2 ], sec^ 1 

p number density of lasing photons, 

p + + p_ , photons ■ cm -3 

p ± number density of lasing photons moving 

along Z-axis, photons • cm -3 

p number density of pumping photons, 

P _ 3 

photons • cm 

p^ normalized pumping photon density 

<J cross section for laser-beam absorption and 

stimulated emission, cm 2 

\\f angle, rad 

[ ] number density, particles • cm -3 

( ) average over cross section of laser tube 

Subscript: 

0 value at Z = 0 

Model Description 


axial distance on which inversion density 
[I*] - [I] /2 > 0 , positive gain length, cm 

active length of laser tube, 33 cm 

length of pump, 1 5 cm 

normalized axial distance downstream of 
Z 

pump entrance, — 

Z L 

Z C+ 

normalized positive gain length, 

Z L 

Z 

normalized pump length, — - = 0.45 

Z L 

constant for C v (eq. (26)), J moF 1 * K -1 
constant for C y (eq. (26)), K 1 
distance along light ray, m 

absorption length for pumping photons, m 

_ 3 

density of lasant gas, kg • m 
angle, rad 

frequency of laser beam, Hz 


Laser Geometry 

The CW laser used in references 14 and 15 is shown 
schematically in figure 1 . The Z-axis of the 1 -D mathe- 
matical model is parallel to the optical axis of the laser 
cavity and points in the direction of the lasant flow. The 
origin is located where the lasant enters the elliptical 
pump chamber. Upstream of this point, the lasant is 
undissociated and does not interact with the laser beam, 
provided the lasant is free of absorbing or scattering 
impurities. The pumping region spans the distance 
0 < Z < Z p , where Z p = 15 cm. We assume that the inci- 
dent pumping radiation at the laser tube is axisymmetric. 
Measurements of the actual incident pumping radiation 
are given in appendix D. Downstream of the pump, the 
lasant has a nonzero inversion density and continues to 
interact with the laser beam until the lasant is withdrawn 
at the end of the tube Z L , where Z-Z L - 33 cm. In the 
computation, the active length of the laser tube 
(0<Z<Z^) is normalized to unity (0 < z < 1), and the 
pump spans the normalized distance 0 < z ^ z p , where 
* p = 0.45. 

Photochemical Reactions 

The kinetic reactions included in the laser model are 
as follows. 


2 



Introduction 

A solar power station advantageously placed in 
space could beam power to other spacecraft and to plane- 
tary surfaces, including the surface of the Earth, as dis- 
cussed in references 1-3. Experiments related to this 
concept include tests of solar-simulator-pumped iodine 
lasers, as discussed in references 4-7. These experiments 
are supported by modeling efforts reported in 
references 8-13. An important finding is that the gaseous 
perfluoroalkyl iodine lasants n-C 3 F 7 I and t-C 4 F 9 I have 
markedly different production rates for the dimers 
(C 3 F 7 ) 2 and (C 4 F 9 ) 2 . Our modeling effort is devoted 
principally to understanding the effects of this differing 
dimer production on laser performance. 

Lee (ref. 14) and Lee et al. (ref. 15) present two 
experimental comparisons of the lasants n-C 3 F 7 I and 
t-C 4 F 9 I when flowed longitudinally in a continuous wave 
(CW) laser. Laser output power is measured in the first 
comparison as the lasant flow speed is varied and in the 
second comparison as the intensity of the solar-simulator 
pump is varied. In both comparisons, the output power 
P oul for t-C 4 F 9 I is found to be about three times greater 
than that for n-C 3 F 7 I. This increase in P oul is not 
explained by an increased utilization of the pump spec- 
trum. For the solar simulator used, the pump spectrum 
utilization is only 20 percent greater for t-C 4 F 9 I than for 
n-C 3 F 7 I. 

The dimer density is represented by [R 2 ], where R 
represents either of the perfluoroalkyl radicals n-C 3 F 7 or 
t-C 4 F 9 and brackets denote the number density. Ershov, 
Zalesskii, and Sokolov (ref. 16) have shown experimen- 
tally that the rate of [R 2 ] production is much less for 
t-C 4 F 9 I than for n-C 3 F 7 I, although a numerical value for 
the ratio of these production rates is not given. Lee 
(ref. 14) notes that the reduced [R 2 ] production increases 
t-C 4 F 9 I recyclability. He also speculates that the reduced 
[R 2 ] production for t-C 4 F 9 I would make more monomer 
radicals [R] available to combine with iodine atoms [I]. 
Consequently, more of the laser lower level [I] would be 
depleted and less [I 2 ], the principal quencher of excited 
iodine [I*], would be formed. This reduction of [I 2 ] 
would also reduce the lasant flow speed. 

Our purpose is to examine these speculations and 
especially to determine the kinetic effects of reduced 
[R 2 ] production on P out . Our approach is to use a one- 
dimensional (1-D) numerical rate equation model. This 
model is described in the next section. The model is 
tuned by using the data given in references 14 and 15; a 
discussion of that process follows the model description. 
Thereafter, general properties of the tuned model and its 
solutions are given. The solutions indicate that the exper- 
imental power curves obtained in references 14 and 15 
could be improved by optimizing the lasant flow speed w 


so that molecular iodine does not build up within the 
pump region of the laser. This optimization is based on 
parametric and scaling studies that are presented in 
appendixes A-D. The flow speed optimization is mod- 
eled in the section “Diagnostic Plots,” and the n-C 3 F 7 I 
and t-C 4 F 9 I lasants are again compared. A fictitious las- 
ant (identical to n-C 3 F 7 I except for a reduced [R 2 ] pro- 
duction rate) is modeled in the section, “Optimization of 
Laser Performance.” The purpose is to isolate the kinetic 
effects of reduced [R 2 ] production. Concluding remarks 
are given next and the appendixes follow. 

Symbols 

B brightness of pump lamp image, 

W -nT 2 rad -1 

C carbon atom 

Cj, C 2 constants of integration 

C specific heat at constant pressure, 

p j-mor'r 1 

c* Cp J kg"' K"' 

MxlO ” 3 

C specific heat at constant volume, 

V r 1 1 

J-moP • K” 

c speed of light in vacuum 

c. three-body, reaction-rate 

l 6—1 

coefficients, cm sec 

_3 

e energy density of pump radiation, J * m 

F incident flux of pump radiation, W/m 

F fluorine atom 

f-C 3 F ? I fictitious lasant identical to n-C 3 F 7 I , except 
that k 3 and k 4 are reduced by factor of 0. 1 

h Planck’s constant, 6.626 x 10 “ 34 J • sec 

I iodine atom 

I p intensity of pump radiation, SC 

I* iodine atom in excited state 

I 2 molecular iodine 

k. two-body reaction rate 

coefficients, cm 3 • sec” 1 

M molecular weight, g mol -1 

P ou{ laser output power, W 

p lasant pressure, Pa 

p ouV laser output power density, w cm -2 



[I*] ) = 4, [RI] + 0.51 [I 2 ] - it, [R] [I*] 

~q 2 { I 2 ] [l*]-ccp( [I*] - | [I] ) 
-9 3 [R] [I*]-? 4 [R 2 ] [I*] 

-? 5 [I*] f I] — *7 [RI] [I*] 

- Cj [RI] [I*] ri] - [RI] [I*] (6e) 

^(w[I]) = 1.49 § 2 [I 2 ] +<?,[ RI1 [I*] 

+ « 2 [I 2 ] [I*] - 2c 5 fl] 2 [R 2 ] 

-* 8 m [Ri] + [Ri] 

+ cap( m - i [I] ) 

- c, [RI] [I*] [I]-2 c 2 [RIJ fl] 2 
+ 2 c 4 [I 2 ] [I] 2 - * 2 [R] [I] 

+ Jfc 4 |RI] [R] + 9 3 [I*] [RI 
+ <? 4 [ I*] [R 2 ] + <7 5 [I*] [I] 

+ k 5 [ R ] [I 2 ] (6f) 

where c is the speed of light in vacuum and p is the num- 
ber density of lasing photons. 

The latter density is given by 

P = P + + P_ (7) 


eliminated from the formulation because equations (8a) 
and (8b) are satisfied by 

p + (Z)p_(Z) = Constant (10) 


The boundary conditions (eqs. (9a) and (9b)) then give 


P _(Z) 


p + (0) 2 

W 7 ) 


and from equation (7) 


p (Z) = p + (Z) + 


P + (0) 2 

R aP+W 


(11a) 


(lib) 


The boundary conditions (eqs. (9a) and (9b)) also reduce 
to 


P + (Z L ) = P + (0) (RR b )~ m (12) 

The formulation is now complete in terms of p + (Z), and 
p_(Z) is determined from equation (11a). The output 
power density p oul in W • cm -2 is given by 

P aa = P + l Z i)V- R J cl " L (13) 

Compressible Fluid Dynamics 

Part of the incident pump power is dissipated by 
heat. The resulting dependence on Z (m) of temperature 
T (K), pressure p (Pa), and speed w (m * sec -1 ) is deter- 
mined approximately from the 1-D, steady-state fluid 
dynamic equations for an inviscid, nonconducting gas 
(ref. 22) as follows. 


where the symbols + and - indicate the direction of pho- 
ton motion along the Z-axis. For steady-state CW opera- 
tion, these photon densities satisfy the equations 

^ = <«a) 

£ --*.<< nVs'i*) < 8b > 

and the boundary condition at Z = 0 

P + (0) = ^P_(0) (9a) 

At Z = Z^ 

p_(Z L ) = R h p + (Z[) (9b) 


Continuity equation . Under all these conditions the 
continuity equation becomes 

= 0 (14) 

where T| is the density (kg ■ m -3 ). Hence, 

C l 

T1 = — (15) 

w 

where the constant C x is given by 

C { = ri 0 w 0 (16) 

where rj 0 is the density and w 0 is the speed at Z = 0. 

Momentum equation. The corresponding momen- 
tum equation is 


where R a and R h are the mirror reflectivities at Z = 0 
andZ=Z L , respectively. The quantity p _(Z) may be 


dw _ 1 dp 

W dZ ~ ~r\dZ 


(17a) 
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Photodissociation reactions . The parent molecule 
RI is irradiated in the pump region, and the following 
reaction occurs: 


RI + hv 

P 




(la) 


where h is Planck’s constant, hv p is the energy of 
pumping photons, is the photodissociation rate of RI, 
and I* = 2 P 1/2 , which is the laser upper level. A sequence 
of reactions produces I 2 , which can also be photodissoci- 
ated in the pump according to the reactions 


^2 

I, + /iv -» 

2 p 

2 

where the photodissociation rate of I 2 , and I = P 3 / 2 , 
which is the laser lower level. Photodissociation rates for 
and £,2 are listed in table I with other laser rate coeffi- 
cients. Data in this table were obtained from 
references 17-21. 


J I + I* (5 1 -percent probability) 
{ 21 (49-percent probability) 


Absorption and stimulated emission reactions . Pho- 
tons in the laser beam can be absorbed by I atoms or 
undergo stimulated emission by excited I atoms accord- 
ing to the reactions 

I* + h v L I + 2 hv L (lc) 

where hv L is the energy of the lasing photons and a is the 
cross section for absorption and stimulated emission. 


Two-body reactions . 

R + I* 4 RI 
*2 

r + i4ri 

R + R -i R, 

* 4 

R + RI -> R, + 1 " 

R + 1 2 4 RI + I 
* 7 

RI + I* — > I 2 + R 

k g 

1 + RI 4 I 2 + R J 

where k, are two-body reaction-rate coefficients. 


( 2 ) 


Three-body reactions. 


I* + 1 + RI 4 I 2 + RI 

C 2 

I + I + RI -> I 2 + RI 

c 4 

I + I + I 2 — ) I 2 + I 2 
c 5 

I + 1 + r 2 — ^ i 2 + r 2 


( 4 ) 


where c t are the three-body reaction-rate coefficients. 

Quenching reactions. The following reactions 
quench the excited state of iodine atoms: 

I* + RI -l I + RI 


I* + I 2 ^I + I 2 

» ?3 

I +R->I + R 

. «4 

I + R 2 ) I + R 2 

* 

I + I-*I + I 
where q ( are the quenching reaction-rate coefficients. 


( 5 ) 


Rate Equations 

For steady-state CW operation, the 1-D rate equa- 
tions are purely functions of Z as shown: 

-|(w[RI]) = *, [R] [I*] +k 2 [ R] [I] 

+ * 5 [R] fl 2 ] -Jfc 7 [I*] [RI] 

-Jk 4 [R] [Rll-5, [ ri 1 

-Jt 8 [I] [RI]-* 9 [RI1 ( 6 a) 

^(w[R]) = ^[RI]-*,rRl [I*] - * 2 [R] [I] 

- 2*3 [R] 2 - [RI] [R] 

- * 5 [R] [I 2 ] + * 7 [RI1 [1*1 
+ * 8 [I] [RI1 + * 9 [RI] 

+ 2i 10 [R 2 ] ( 6 b) 

-^(w[R 2 1) = ^ 3 [R] 2 + A' 4 [RI] [R] -*, 0 [R 2 ]( 6 c) 


Pyrolysis . The lasant RI and dimer R 2 can also 
undergo thermal dissociation, especially at high tempera- 
tures. Thus, 


RI -4 R + I 
*10 

r 2 ->r + r 


( 3 ) 


^(w[I 2 ]) - c, [RI] [I*] [I] +C 2 [RI] [I ] 2 
+ c 4 [I 2 ] [I] 2 -^ 2 [I 2 ] 

+ Jt 7 [RI] [I*] ~k 5 [R] n 2 ] 

+ c 5 [I] 2 [R 2 ] + ^ 8 [RI] [I] ( 6 d) 
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experimental conditions with these most favorable rate 
coefficients). Therefore, some adjustment of the rate 
coefficients was necessary for the theoretical study. Most 
rate coefficients are not measured directly but are 
inferred semiempirically through a combination of mea- 
surements and models. Thus, uncertainties in published 
values of rate coefficients depend not only on measure- 
ment uncertainties but also on model uncertainties. The 
latter uncertainties often are not evaluated. 

Wishing to depart from the published values as little 
as possible, we performed a sensitivity study with the 
model. This study showed that under the given experi- 
mental conditions, the model was most sensitive to k 2 , £ 3 , 
& 5 , and q 2 and the output power was increased by 
increasing k 2 and k$ and by decreasing k 3 and q 2 . 
Accordingly, we multiplied k 2 and k 5 and also the pump 
photodissociation rates ^ and £ 2 by a factor of 2.55 and 
divided and q 2 by the same factor to achieve lasing 
and obtain the n \ -model curve of figure 2(a). These 
adjusted values for n-C 3 F 7 I are listed in table I, where 
temperature dependence is included where known. 

The /I -model curve of figure 2(a) was obtained next 
by using the t-C 4 F 9 I data given in table 11(a) and the cor- 
responding rate coefficients given in table I. As shown in 
table I, the principal difference assumed for t-C 4 F 9 I is 
that the production rate of the dimer fR 2 ] is reduced by a 
factor of 0.1 (reactions k 2 and & 4 ). The other rate coeffi- 
cents and parameters are taken to be the same as those for 
n-C 3 F 7 T, except for M, a, (3, and Various factors for 
reduced [R 2 ] production were also tried, but the factor 
0.1 gave the best match, subject to the values assumed 
for the other rate coefficients of t-C 4 F 9 I. 

We had originally found that the n2~ and /2-model 
curves almost overlapped. Recalling that chemical analy- 
sis of the n-lasant had revealed a 20 -percent dimer impu- 
rity, as shown in appendix A, we incorporated this 
feature into the model. This change succeeded in increas- 
ing somewhat the vertical separation of the theoretical 
curves in figure 2 (b) without much affecting those in fig- 
ure 2(a). Chemical analysis of the t-C 4 F 9 I lasant (appen- 
dix A) indicated an even higher amount of an unknown 
impurity; however, a 20 -percent dimer contaminant for 
both n-C 3 F 7 I and t-C 4 F 9 I gave the best match in figure 2. 
This adjustment completed our tuning of the model. 

Diagnostic Plots 

Besides P oul , the laser model computes and plots the 
following densities as functions of c: the parent molecule 
[RI], the monomer [R], the dimer [R 2 ], molecular iodine 
P 2 ], excited iodine [I ], ground-state iodine [I], the inver- 
sion density [I*] - [I]/2, and the lasing photon densities 
p + and p_. Also computed and plotted are T , f|, p , and w. 
The plots of these computed quantities (diagnostic plots) 


are useful in the analysis of laser performance. 
Figures 3-6 present diagnostic plots for the four model 
curves of figure 2 . 

Diagnostic plots for both end points of the n \ -model 
curve of figure 2(a) are shown in figure 3. These diag- 
nostic plots show that the laser is not operating opti- 
mally. In figure 3(a), [I 2 ] rapidly builds up at the 
midpoint of the pump region (positive normalized gain 
length £ C+ ) and remains high downstream of this point 
until, at the exit of the pump region ( z p ), it reaches an 
even higher plateau. The effects of this [I 2 ] buildup on 
[I ], [I], and the inversion density [I ] - [I]/2 are shown 
in figures 3(b)-3(d). In particular, the inversion density 
becomes negative midway through the pump region and 
remains so downstream of this point, although it 
becomes somewhat less negative at the exit of the pump 
region. 

These modeling results are consistent with the well- 
known observations (ref. 19) that [I 2 ] is a strong 
quencher of [I ] and can make the inversion density 
become negative. This negative inversion density results 
in a net absorption of lasing photons on the downstream 
half of the pump section, as shown by the curves for the 
lasing photon densities p + and p_ in figure 3(e). In partic- 
ular, p + reaches its peak value in the middle of the pump 
section and thereafter decreases all the way to the output 
mirror at z = 1. Because p oul is proportional to p + at the 
output mirror, the model laser is not operating optimally. 
For optimal operation, p + would continue to increase 
until the flow exits the pump at z = 0.45. 

Other features of this model run warrant explanation 
as well. Because [R] reacts much more readily with fl] 
than with [I ] to form [RI] (i.e., k 2 » &j), [R] decreases 
rapidly as [I] increases at the midpoint of the pump sec- 
tion, as shown in figure 3(f). The distribution of [R 2 ] is 
nonzero at z = 0 because of the 20 -percent initial [R 2 ] 
contaminant, as shown in figure 3(g). 

Figures 3(h)-3(k) show plots of T, w, rj, and p 
(almost constant). These four plots are representative of 
all the model runs and are subsequently omitted. The 
moderate increases in T and w from solar-simulator- 
pump heating generally agree with the experiments in 
which a water-cooled quartz laser tube was used. (See 
fig. 1 .) 

The t \ -model curves of figure 4 and one /i2-model 
curve {I p =1100 SC) of figure 5 also show similar non- 
optimal laser performance. However, in the other 
«2-model curve ( I p = 450 SC) of figure 5 and in the 
/ 2 -model curves of figure 6 , the buildup of [I 2 ] occurs 
only downstream of the pump, which indicates optimal 
operation. 
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From equation ( 1 6), this equation may be written as 


rfZ <p + C = 0 

(17b) 

Upon integration, we obtain 


p = C 2 - Cj w 

(18) 

where 


^2 “ ^0 ^ 

(19) 


heat at constant volume C v (J • mol 1 ■ K ! ) is closely 
approximated by 

C v = a exp [p (T- 300)] (26) 

where a and f} are constant for each lasant as given in 
table I. The specific heat at constant pressure C p in 
J ♦ moF 1 * K _1 is 

C = C +R (27) 

p v v ; 

Thus, we have (in J • kg -1 • K -1 ) 


Equation of state . If the lasant is idealized as a per- 
fect gas, the equation of state is given by 


C* = — -exp[3(T-300)]+R* (28) 

Mx 10 


p = 0 R*T 


(20) and from equation (22) 


Here, R* (J • kg 1 • K ! ) is given by 

R* = (21) 

Mx 1CF 3 


where the gas constant R = 8.314 J • mol -1 - K _l , and 
M(g ■ moF r ) is the molecular weight of the lasant. Solv- 
ing equation (20) for T, we obtain 


T 


vv 

7 




( 22 ) 


From equation (20), C x may also be written 


C, = 


Pq w 0 

R*T n 


(23) 


where p 0 , w o> an< ^ 310 measured quantities. 


Energy equation. The 1 -D, steady-state energy 
equation for an inviscid, nonconducting perfect gas is 
given by 


dT dp 

c p r ' w Tz- w dz m 2 


(24) 


where C* (J ■ kg -1 * K _1 ) is the specific heat at constant 
pressure divided by Mx 1CT 3 , and Q (W * m -3 ) is the 
heat from the incident pump radiation. From 
equations (15), (18), and (22), this equation may be 
rewritten as 


dw _ R*Q 

dZ R*C l w + C*(C 2 -2C ] w) 


(25) 


Specific heats. In equation (25), C* may be approxi- 
mated from data in reference 22 as follows. The specific 


C = 

u p 


a 


Mx 10 


-3 


exp 


r 

r 

f r ^ 

- 

1 

p 

W 

2 

— — vv 

-300 


l 

,R* 

l C l J 


1 


' + R (29) 


Equation (29), which for a given lasant gas gives C* as a 
function of vv, is appropriate for substituting back into 
equation (25). 


Numerical integration procedure . For a given heat- 
ing rate Q(Z) and values for the constants Cj and C 2 
(obtained by measuring p 0 , T 0 , and w 0 at Z = 0), 
equation (25) can be numerically integrated to give vv(Z). 
The other fluid dynamic fields, r\ (Z), p(Z), and T(Z) fol- 
low from equations (15), (18), and (22), respectively. 
The flow speed w(Z) appears explicitly in the rate 
equations (6a)-(6f). The density r\(Z) determines the par- 
ent molecule number density [RI], and T(Z) enters the 
temperature-dependent rate coefficients. 


Model I\ining 

The rate equation model was tuned by matching as 
closely as possible the experimental curves for the las- 
ants n-C 3 F 7 I and t-C 4 F 9 I, as given in references 14 and 
15. The laboratory data for these curves are given in 
tables 11(a) and 11(b), where the original lasant flow units 
in standard cubic centimeters per second have been con- 
verted to meters per second. The wide range of pressure 
and pump l p intensity values makes these data sets 
appropriate for tuning the model. The best match 
achieved with the model is shown in figure 2. This 
match, although imperfect, was obtained with difficulty 
as described in the discussion that follows. 

The data in table 11(a) for the n 1 experiment curve of 
figure 2(a) gave a theoretical P out = 0 W when used with 
the published rate coefficients for n-C 3 F 7 I that are the 
most favorable for lasing (i.e., the model did not even 
reach the threshold for lasing under the given 
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reduced [R 2 1 production are highly beneficial to efficient 
operation of a CW iodine laser. 

Concluding Remarks 

Two experiments which compare the performance of 
the lasants n-C 3 F 7 I and t-C 4 F 9 I can be approximated by a 
one-dimensional numerical rate equation model. In this 
model, the principal difference assumed for the lasants is 
that the dimer production rate for t-C 4 F 9 I is one order of 
magnitude less than for n-C 3 F 7 I. The model results indi- 
cated that laser output power could be increased by opti- 
mizing the lasant flow speeds so that the principal 
quencher, molecular iodine, does not build up within the 
pump region. This optimization method was based on 
parametric and scaling studies that are also presented. 
Such optimized model runs showed that t-C 4 F 9 I had a 
larger output power than n-C 3 F 7 I, although the power 
increment depended on the operating conditions. 

Optimization of the flow speeds improved the basis 
for comparing n-C 3 F 7 I and t-C 4 F 9 I; however, the kinetic 
effects of reduced dimer production were still not clear 
because the two lasants had different molecular weights, 
pump spectrum utilizations, densities, and flow speed. To 
clarify these kinetic effects and determine how reduced 
dimer production results in greater output power, we also 


modeled a fictitious lasant f-C 3 F 7 I, which differed from 
n-C 3 F 7 I only by a reduced dimer production rate. The 
results of this theoretical comparison confirmed earlier 
speculations; that is, simply by reducing dimer produc- 
tion, we produce more monomer radicals (C 3 F 7 ). These 
radicals then combine with iodine atoms to enhance 
depletion of the laser lower level and reduce the growth 
of the principal quencher, molecular iodine. Both of 
these effects tend to increase the lasing photon density 
and, hence, the output power. 

This theoretical study also found that an order-of- 
magnitude decrease in the dimer production rate halves 
the lasant flow speed required to prevent the buildup of 
molecular iodine in the pump region. Also, more lasant 
molecules are recovered after they pass through the laser 
tube and less iodine molecules are produced; hence, 
reduced dimer production improves lasant recyclability. 
If other properties are equal, reduced dimer production is 
clearly a desirable feature for the flowing lasant in a con- 
tinuous wave iodine laser. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
December 23, 1994 
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Optimization of Laser Performance 

Although we cannot prevent the buildup of [I 2 ] 
downstream of the pump, we can prevent it inside the 
pump region, where it does the most harm, by increasing 
the lasant flow speed w. This effect of vv — together with 
other factors that could be used to control the buildup of 
[I 2 ] — is discussed in appendix B. In this context, we take 
the optimal w to be the minimum flow speed for which 
the buildup of [I 2 ] occurs only downstream of the pump. 

We found that optimizing the theoretical power 
curves in figure 2 provides a better way to compare the 
lasants n-C 3 F 7 I and t-C 4 F 9 I. Optimizing the n i- and 
t\ -model curves of figure 2(a) by increasing w and by 
eliminating the initial dimer contaminant (i.e., setting 
[R2(0)] = 0) gives the increased P out shown in figure 7. 
These curves show that under optimal conditions, 
t-C 4 F 9 I gives about twice as much output power as 
n-C 3 F 7 I at about half the flow speed. Figures 8 and 9 
present diagnostic plots for both end points of each opti- 
mized curve shown in figure 7. These plots confirm that 
the buildup of [I 2 ] occurs downstream of the pump and 
that the inversion density fl ] - [11/2 is positive through- 
out the pump region. 

Similarly, optimization of the n2- and /2-model 
curves of figure 2(b) gives the increased P out shown in 
figure 10. Again, t-C 4 F 9 I has higher output power than 
n-C 3 F 7 I, although the results are less dramatic. We con- 
clude, for optimized flows, that t-C 4 F 9 I is superior to 
n-C 3 F 7 I; however, the difference in output power 
depends on the operating conditions. 

Kinetic Analysis of Reduced Dimer Production 

Although optimization of the power curves, as 
shown in figures 7 and 10, facilitates comparison of the 
lasants n-C 3 F 7 I and t-C 4 F 9 I, it does not clarify the kinetic 
effects of reduced dimer production because the opti- 
mized /-model curves — besides having reduced dimer 
production relative to the optimized /i-model curves — 
also have different pump-spectrum utilizations, molecu- 
lar weights, specific heats, densities, and flow speeds. 

To analyze unambiguously the effects of reduced 
[R 2 ] production, we considered a fictitious lasant f-C 3 F 7 I 
in which the only difference from n-C 3 F 7 I is reduced val- 
ues for k 3 and & 4 , We start with the n \ -model optimal 
curve of figure 7. Because [I 2 ] builds up only down- 
stream of the pump, the same will be true for the 
f\ -model (obtained from the n \ -model optimal curve by 
reducing k 3 and k 4 by a factor of 0. 1 without changing vv 
or anything else). 

Reduced dimer production (without other changes) 
gives the increased ^out shown by the f\ -model curve 


relative to the n \ -model optimal curve in figure 1 1 . Diag- 
nostic plots for the /I -model are presented in figure 12 
for comparison with the diagnostic plots for the 
n \ -model optimum, as presented in figure 8. These diag- 
nostic plots are supplemented by detailed plots of [I 2 ] 
and [RI] for the two models in figure 13. 

These diagnostic and detailed plots show clearly that 
a reduced [R 2 ] production rate decreases [R 2 ], [I 2 ], [I], 
[I ], and [I ] - [Tj/2 in the pump region and increases [R], 
[RI], p + , and p_. The results confirm the speculations 
made in reference 1; that is, the decrease in [R 2 ] and [I 2 ] 
and the increase in [RI] improve lasant recyclability. The 
reduced [R 2 ] production results in a greater density of 
monomer radicals [R] in the pump region that can com- 
bine with iodine atoms [I] and thus enhance depletion of 
the laser lower level and reduce the buildup of the princi- 
pal quencher [I 2 ]. The lower value of [I 2 ] increases p + 
and p_ because fewer downward transitions occur by 
quenching and more by the stimulated emission of lasing 
photons. The lower value of [I] reduces the absorption of 
p + and p_. Hence, the reductions in [I 2 ] and [I] both help 
increase the output power P oul , which is proportional to 
p + at the output mirror. 

After repeating this reduced dimer analysis with the 
n2-model optimal curve of figure 10 and the fictitious 
lasant f-C 3 F 7 I, we obtain the increased P out shown by the 
fictitious /2-modcl curve in figure 14. Figures 1 1 and 14 
clearly show that simply by reducing the dimer produc- 
tion rate, we were able to increase output power, 
although the amount of increase depends on the operat- 
ing conditions. 

Although not shown, the diagnostic plots for the 
^2-model optimum and the /2-model are similar to those 
given in figures 8 and 12 for the n \ -model optimum and 
the f\ -model, respectively. The corresponding detailed 
plots are given in figure 15, which again shows that the 
reduction of [R 2 ] production increases [RI] and decreases 
[I 2 ]. This second comparison reinforces the previous 
conclusions about the benefits of reduced [R 2 ] 
production. 

Finally, we investigate the speculation that reduced 
[R 2 ] production allows the lasant flow speed w to be 
reduced. In figure 11, the values of vv 0 for the fictitious 
/I -model are the same as for the n\ -model optimum. We 
now optimize the f\ -model by reducing vv 0 to the mini- 
mum values that still prevent the buildup of [I 2 ] inside 
the pump. We find that Wq can be halved and that P out is 
not significantly affected, as shown in figure 16. Diag- 
nostic plots for the f\ -model optimum (halved vvq) are 
shown in figure 17. These plots confirm that [I 2 ] builds 
up only downstream of the pump. Although not shown, 
similar results are obtained when the fictitious /2-model 
of figure 14 is optimized by halving vv 0 . These effects of 
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Appendix A 

Gas Chromatographic-Mass Spectrometric 
Analysis of Lasants n-C3F 7 I and t-C 4 F 9 I 

A gas chromatographic-mass spectrometric (GCMS) 
device was used to examine samples of the two lasants. 
The findings are summarized below. 

Analysis of n-C3F 7 I Sample 

No significant difference was noted between the pre- 
and postlased n-C 3 F 7 I material. The analysis revealed the 
following compounds: 

(C 3 F 7 ) 2 -20 percent 

C 2 F 5 I Trace amounts 

n-C 3 F 7 I =80 percent 

Unknown perfluorocarbon Trace amounts 

C 2 F 4 C1I Trace amounts 


The presence of I 2 was not detected in the GCMS 
analysis, even though I 2 was visually detectable in the 
postlased sample. Several explanations can be given for 
the absence of I 2 in the postlased sample: the I 2 peak was 
masked by either the (C 3 F 7 ) 2 dimer peak or the n-C 3 F 7 I 
peak, the I 2 was trapped on the gas chromatograph col- 
umn, or the I 2 concentration was so small that it was not 
detected on the gas chromatograph. 

Analysis of t-C 4 F 9 I Sample 

More than 20 gas chromatographic peaks were 
detected. Only the ethyl alcohol and t-C 4 F 9 I peaks could 
be identified. The remaining peaks contained perfluoro- 
alkyl iodine compounds and perfluoroalkanes. The major 
peaks are listed below and are approximate. 


Peak 

Identification 

Sample, percent 

15 

Ethyl alcohol 

3 

62-83 

t-C 4 F 9 I 

35 

444-480 

Unknown 

60 

1154 

Unknown 

2 
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Appendix B 

Parametric and Scaling Studies 

As shown in figure 3, the inversion density [I*] - [I]/ 
2 > 0 for only a fraction of the pump length Z p . We 
termed this fractional distance the positive gain length 
Z G+ - We performed a parametric study to determine how 
Z G+ varies with flow speed w, number density [RI], and 
pump intensity l p . The result is expressed by the 
proportionality 


The output power ^out in W scales with the inner 
radius r t of a circular cylindrical laser tube according to 
the proportionality 

P ou|-<IVAt>'-f < B3 > 

where { p p N ) is the normalized density of pumping pho- 
tons averaged over the cross section. A derivation of 
(P p n) as a function of the ratio r/ 8, where 5 is the pump- 
ing photon absorption length, is given in appendix C. 



An increase in w retards the buildup of [I 2 ], whereas 
an increase in [RI] and I p increases the level of 
photoionization and accelerates the buildup of [I 2 ]. This 
proportionality was used to optimize the laser model (to 
achieve positive gain throughout the pump region) by 
increasing w until z G+ = z p . 

A scaling study was also performed to determine 
how output power density p out (W • cm -2 ) scales for sim- 
ilar lasers — lasers that have the same ratio of positive 
gain length to pump length Z G JZ p . For such similar 
lasers, p out was found to scale according to the 
proportionality 

/\> ut ~ z G+ y R n < B2 > 


The scaling study may be summarized as follows: 
Two optimal CW iodine lasers are similar if the lasant 
flow speeds satisfy the relation 


^2 _ Z p2 2-Fp2 

w i ~ z P i [RI] iVti 


(B4) 


For such similar lasers, the output power P out scales as 

P o»u2 _ Z p2 t Rl 2 j>2 r j2 (PpNl) 

P «,l ~ Z „ ™ rf, <P^1> 

A plot of (p pN ) as a function of rj 5 is given in 
appendix C. 
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Appendix C 

Derivation of Average Pumping Photon 
Density in Absorbing Lasant With Circular 
Cylindrical Symmetry 

This appendix provides a derivation of the quantity 
(P p n) used to scale the output powers of equation (B5) in 
appendix B. 

Two-Dimensional Relation of Incident Flux to 

Transmitting Surface Brightness 

Assume that the laser tube and the incident pump 
radiation are axially symmetric and independent of Z. 
(Appendix D gives a discussion of how this idealization 
compares with theoretical and actual elliptical pump 
chambers.) We also assume that the optical image of the 
pump source fills the interior of the laser tube, as shown 
in figure 18, where the curved wall of the circular glass 
tube acts as a negative lens. Furthermore, we assume that 
in the absence of absorption the pump intensity is uni- 
form and isotropic on a cross section of the pump image. 
Then, if we take a point on the perimeter, as shown in 
figure 19, the image brightness B (W ♦ m“ 2 • rad" 1 , inter- 
preted as watts per meter perimeter per meter depth per 
radian) can be integrated to determine the incident flux of 
pump radiation F (W • m~ 2 ): 

r /2 

B cos 0 dB = 2 B (Cl) 

71/ 2 

where the angle 0 is in radians. A comparable 3-D rela- 
tion is given in reference 23 on page 23. 

Pumping Photon Density Transmitted by Surface 
Element 

We now determine the differential pumping energy 
density de (J • m -3 ) at an interior point at depth y due to 
incident radiation passing through length dS of the 
perimeter, as shown in figure 20. Initially, assume that 
the laser tube is evacuated so that absorption is not a fac- 
tor. Then, the power within ray dB (W • nT 1 depth) 
is given by B cos 0 dS dB. Because this power results 
from photons moving at the speed of light in vacuum c , 
the pumping energy (J • m“ 2 (per meter depth per meter 
length)) within dB is given by ( B/c ) cos 0 dS dB. But 


dB = dsly ; hence, de (J • m 3 ) due to the incident radia- 
tion that penetrates dS is given by 

de = — cos0 dS (C2) 

cy 

and the corresponding pumping photon density dp p 
(photons • itT 3 ) is given by 

dp = — — cos9 dS (C3) 

r P hv cy 
P 1 

For an absorbing lasant gas, this photon density becomes 
dPp ~ hv c cos0 dS (C4) 

VpCj 

where 5 is the absorption length at the pump frequency. 

Total Pumping Photon Density at Interior Point 
Due to Axisymmetric Incident Flux 

Substitution of equation (Cl) into (C4) gives de at 
point H of figure 21 when F is transmitted through dS : 

dpp = 2 j ~ — exp(-y/5) cos0 dS (C5) 
nv p cy 

However, by introducing the angle \\f shown in figure 21 
we may write 

dy = dS (C6) 

and equation (C5) becomes 

dp p = exp ( ~ y/5) d v (C7) 

n\p 

_ T 

The total pumping photon density (photons • m ) at 
point H in an absorbing gas may then be written 

= j“e^p ( -Y/S ) rf V (C8) 

p u 

and by the law of cosines 

[2 2 Y~ 

y = -rcos\|/+ Ar t -r sin \| i (C9) 
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Hence, the pumping photon density (photons • m 3 ) at 
radius r in an absorbing lasant is given by 


f’p lr - r r i > - 7^7 f>P 


r< 

J o 


4 


2 2 2 

r f -r sin y 


1 

LB 


dy 


/■cosy 


(CIO) 


Using this quantity to normalize the photon density in 
equation (CIO), we obtain 


P pN^ r N 


-v 5) = i r 

' Jt Jo 


2 2 
V sin y 


eX P|_5 

d y 


(Cl 3) 


The average pumping photon density over the cross sec- 
tion is given by 

(P p ( r r 5 >> = f'dr2nr p p (r,r r 5) (Cl 1) 
Kr" { 0 

Normalized Pumping Photon Densities 

In the limit as 5 — » °o s equation (CIO) gives the uni- 
form photon density for a nonabsorbing lasant as 

p c (5 ^” )= ^f <C12) 

p 


where = rlr t is the normalized radius. Plots of this nor- 
malized pumping photon density versus the normalized 
radius are given in figure 22 for various values of r / 5 . 

The average pumping photon density in an absorbing 
lasant, as given by equation (Cl 1), can also be normal- 
ized by 

(P pN (r/S)) = 2^dr N r N p pN {r N ,r t !?>) (Cl 4) 

This normalized average pumping photon density is plot- 
ted versus r I 5 in figure 23. This quantity is an important 
factor in the scaling relation (eq. (B5) from appendix B) 
for calculating how output power scales with the radius 
of the laser tube. 
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Appendix D 

Theoretical and Experimental Performance of 
Elliptical Pump Chamber 

A ray tracing for the elliptical cross section of the 
pump chamber used in the laser experiments is shown in 
figure 24(a). The rays are assumed to emanate isotopi- 
cally from a line source (perpendicular to the page) at the 
right focus. However, the rays incident on the laser tube 
at the left focus are found to be nonisotropic. As shown 
in figure 24(a), they are more concentrated on the side 


toward the source. A decrease in the distance between 
the focii in the elliptical pump chamber makes the inci- 
dent radiation at the laser tube more isotropic, as shown 
in figure 24(b). 

An experimental polar plot of the incident pumping 
radiation measured inside the laser tube is shown in fig- 
ure 25. Differences from the theoretical plots result from 
the finite size (1 1 mm diameter) of the pump lamp, from 
obstruction by the lamp start-up wire, and from imper- 
fections in the curvature and finish of the optical sur- 
faces, including the elliptical cylinder and the flat end 
plates. 
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Table I. Lasant Rate Coefficients and Other Parameters 




Lasant 



Parameter 

Unit 

n-C^F^I 

t-C4F9l 

Reference 


3 -1 

cm • sec 

1 x 10" 14 

Same 

17 

% 

3 -1 

cm ■ sec 

(2.04 ± 1.02) x 10 -11 

Same 

18 

% 

3 -1 

cm • sec 

(7.84 ± 1.18) x 10" 13 

(7.84 ± 1.18) x 10“ 14 

18 

k 4 

3 -1 

cm • sec 

3 x 10" 16 

3 x 10 -17 

19 

% 

3 -1 

cm • sec 

2.55 x 10“" 

Same 

19 

k 7 

3 -1 

cm • sec 

(3 ± 1.5) x 10" 19 

Same 

19 

*8 

3 -1 

cm • sec 

1.6 x 10" 23 

Same 

19 

k 9 

cm 3 • sec -1 

1 x 10 15 exp(-2.48 x I0 4 /T) 

Same 

19 

*10 

3 -1 

cm ■ sec 

1 x 10 17 exp(-4.73 x \0 4 /T) 

Same 

19 

Q\ 

3 -1 

cm * sec 

(1.7 ± 0.2) x 10 -17 

Same 

18 

*92 

3 -I 

cm • sec 

(1 .49 ± 0.10) x 10 _1 1 

x exp [-4. 4 x 1 0 _3 (7-300)] 

Same 

19 

93 

3 -1 

cm ■ sec 

3.7 x 10~ 18 

Same 

19 

94 

3 -I 

cm ■ sec 

4.7 x 10“ 16 

Same 

19 

95 

3 -1 

cm • sec 

1.6 x 10 -14 

Same 

19 

C\ 

cm 6 ■ sec” 1 

1.6 x 10 -33 

Same 

20 

c . 2 

cm 6 ■ sec” 1 

(5.7 ± l)x 10~ 33 exp[(1360±200)/71 

Same 

19 

c 4 

cm 6 ■ sec -1 

antilog 10 {-29.437 - 5.844 log 10 (77300) 
+ 2.163 [log 10 (77300)] 2 } 

Same 

19 

c 5 

cm 6 ■ sec” 1 

(8 ± 2) x 1 0~ 33 exp[( 1 3 1 0 ± 1 00)/7] 

Same 

19 

% 

sec” 1 

(1.2 ± 0.4) x 10 _2 / p 

(1.44 ± 0.48) x 10"% 


% 

sec” 1 

(1.2 ± 0.4) x 10 _! / p 

Same 


a 

J ■ mol -1 • K” 1 

147.23 

183.26 

21 

P 

K- ] 

1.2 x 10' 3 

1.4 x 10" 3 

21 

M 

g • mol” 1 

296 

346 


[R 2 (0)] 

-3 

cm 

0.2 [RI (0)] 

Same 



*Cited value multiplied by 2.55. 
t Cited value divided by 2.55. 
^Semiempirical value multiplied by 2.55. 
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Figure 1. Laboratory schematic of CW iodine laser with longitudinally flowing lasant gas and continuous pumping by 
argon arc lamp. 
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(a) Data from reference 14. 



(b) Data from reference 1 5. 


Figure 2. Experimental comparison of output power from n-C 3 F 7 l and 


t-C 4 F 9 I with corresponding model-tuning curves. 
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(k) Lasant pressure p versus z. 
Figure 3. Concluded, 
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(g) Perfluoroalkyl dimer density [R 2 ] versus 
Figure 5. Concluded. 



Pump region 


z 


(a) Molecular iodine density [y versus z. 



(b) Excited iodine density [I*] versus z. 

Figure 6. Theoretical plots of diagnostic quantities versus z for both end points of /2-model curve of figure 2(b). 





(c) Ground-state iodine density [I] versus z. 



Pump region z 

(d) Inversion density [I*] - [I]/2 versus z. 
Figure 6. Continued. 
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[R], Particles • cm ° D . . • 

1 p +) Photons • cm 


i 


Pump region z 

(e) Lasing photon densities p + and p_ versus z 



(f) Perfluoroalkyl radical density [R] versus z. 


Figure 6. Continued. 


Pump region 


z 


(g) Perfluoroalkyl dimer density [R 2 ] versus 


Figure 6. Concluded. 





Figure 7. Theoretical effects on output power of optimized lasant flow speeds in n\- and fl -models of figure 2(a). 
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(a) Molecular iodine density [I 2 ] versus z. 



(b) Excited iodine density [I*] versus z- 

Figure 8. Theoretical plots of diagnostic quantities versus z for both end points of n \ -model optimal curve of figure 7. 
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(a) Molecular iodine density fl 2 ] versus c. 



(b) Excited iodine density [I*] versus z- 

Figure 9. Theoretical plots of diagnostic quantities versus z for both end points of fl -model optimal curve of figure 7 




(c) Ground-state iodine density [I] versus z. 



Pump region 


(d) Inversion density [1*1 - [I]/2 versus z. 
Figure 9. Continued. 
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(g) Perfluoroalkyl dimer density [R 2 ] versus z . 
Figure 9. Concluded. 
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n2-model 

n2-model optimum 
t2-model 

t2-model optimum 


Figure 10. Theoretical effects on output power of optimized lasant flow speeds in n2 - and r2-models of Figure 2(b). 



Figure 11. Theoretical comparison of output power from fictional/ 1 -model and n \ -model optimum of figure 7. 




(a) Molecular iodine density [I 2 ] versus z. 



(b) Excited iodine density [I*] versus z . 

Figure 1 2. Theoretical plots of diagnostic quantities versus z for both end points of fictional/ 1 -model curve of figure 1 1 . 


47 












[Rl], Particles • cm [IJ, Particles • cm 


nl -model optimum 
fl -model 



(a) Molecular iodine density [I 2 ] versus Z- 



(b) Lasant number density fRI] versus z. 


Figure 13. Detailed theoretical plots of Fictitious/1 -model and n\ -model optimum of figure 1 1 at w 0 = 19.5 m • sec 





Pump region 


(b) Lasant number density [RI1 versus z. 

Figure 15. Detailed theoretical plots for fictitious/2-modeI and ^2-model optimum of figure 14 at I n = 740 SC. 


out’ 



Figure 16. Theoretical comparison of output power from n I -model optimum and fictitious/1 -model of figure 1 1 with 
/I -model optimized by halving the lasant flow speeds. 





(a) Molecular iodine density [LJ versus z* 



(b) Excited iodine density [I*] versus z . 

Figure 17. Theoretical plots of diagnostic quantities versus z for both end points of fictitious/1 -model optimal curve of 
figure 16. 
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*— Pump light image Pump light 

in laser tube 

Elliptical cavity 

Figure 1 8. Cross section of elliptical pump chamber* 



Figure 19. Cross section of laser tube with pump image brightness B and incident flux F. 










tube source 

(a) Elliptical section used in laser experiments. 



Laser tube 


Line source 


(b) Elliptical section with closer focii. 

Figure 24. Theoretical ray tracing for elliptical pump chamber with line source at right focus. 
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Figure 25. Polar plot of normalized incident pumping intensity I p measured at inner radius of laser tube for elliptical 
pump chamber used in laser experiments. Corresponding theoretical ray tracing is shown in figure 24(a). 
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